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Directly interacting particles are considered in the multitime formalism of predictive rela- 
tivistic mechanics. When the equations of motion leave a phase-space volume invariant, it 
turns out that the phase average of any first integral, covariantly defined as a flux across a 
7 n- dimensional surface, is conserved. The Hamiltonian case is discussed, a class of simple 
models is exhibited, and a tentative definition of equilibrium is proposed. 
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1. Introduction 

It is a pleasure to contribute this paper in the honor of prof. Larry P. Horwitz who has 
done (and keeps doing !) so much in practically all the domains of relativistic dynamics. 

Relativistic statistical mechanics has a long story, but we may notice that, whereas 
the theory of relativistic ideal gases has received deep and detailed developments, little 
has been achieved in order to account for mutual interactions beetwen particles. Most 
work available in the litterature concerns kinetic theory [1], and it is only through wild 
simplifications that collision terms can be somehow "derived" from the underlying n-body 
dynamics. For instance various results are obtained with the help of a straight-line ap- 
proximation which amounts to calculate the force applied to each particle as if all other 
particles were undergoing free motion. This treatment retains very little information (if 
any) about the underlying dynamics. 

Another approach consists in giving up manifest covariance in order to deal with a more 
familiar framework, where the concepts and methods of newtonian statistical mechanics 
are intensively employed. This approach is plagued with unpleasant complications (and 
sometimes serious difficulties) concerning relativistic invariance. 

Looking for a more satisfactory theory, where interactions between particles would be 
described in a fully covariant fashion, several authors [2] [3] have hesitated between two 
points of view: field theory and action at a distance. 

In fact, if one aims at rigorous results, it is hardly possible to deal with the tremendous 
infinity of the field degrees of freedom. It seems more convenient to eliminate the field 
variables, and to formulate equations of motion for the particle degrees of freedom only. 
This is the case for Feynman- Wheeler electrodynamics [4] , or more generally for any inter- 
action a la Van Dam and Wigner [5]. But these versions of relativistic dynamics are not 
yet tractable; in the best cases, they lead to difference-differential equations (they involve 
retarded arguments) and seem to require an infinite number of initial data for determining 
the trajectories. 

In contradistinction, the modern theories of directly interacting particles [6] are tailored for 
dealing with only a finite number of degrees of freedom. They naturally suppose that the 
field degrees of freedom have been eliminated, but in addition they assume from the start 
that a finite number of initial data is sufficient in order to determine the evolution of the 
system. In this sense they are a simplification of the approach initiated by the authors of 
refs [4] [5]. 

In this spirit, we have proposed many years ago that an n-body relativistic system be 
described by truly differential equations of motion, involving n independent evolution 
parameters [7] [8]. Our departure from the Feynman- Wheeler scheme is the requirement 
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that the right-hand sides of the equations depend only on the positions (in space-time) and 
their first derivatives [9] . Following the specialized litterature, let us refer to this approach 
as predictive mechanics; clearly this terminology accounts for the possibility to determine, 
in principle, world lines from the knowlege of (a finite number of) initial data [10]. 
Lapiedra et al. [11] have undertaken an application of predictive methods to relativistic 
statistical mechanics. However their definitions for averages, normalization and equilibrium 
are not manifestly covariant. 

A quite different view of direct relativistic interactions has been proposed by Horwitz 
and Piron [12]; in their approach, a unique "historical time" is supposed to describe the 
motion of any number of particles; in other words a single-parameter evolution group acts 
in a finite-dimensional covariant phase space. Along this line, Horwitz et al. [13]have 
extensively developed a "statistics of events" considered off the mass shell. 
In contradistinction we intend to build a statistics of orbits, in agreement with Hakim's 
philosophy (which is after all very close to the principles of predictive mechanics) . 

In a covariant framework like ours, the relativistic analog of a function integrated over 
phase space usually appears as the flux of a generalized current accross some manifold [14] 
[2]. For instance the n-body distribution function N is normalized through the formula 
J Na — 1, where a is, in our notation, a 7?i-differcntial form introduced by Hakim [Na 
is interpreted as the numerical Hux of particles) . Independence of the result with respect 
to the manifold where integration is performed crucially requires that the integrand be a 
closed form in the sense of differential geometry. 

In fact the importance of having d{Na) = was clearly emphasized in ref. [2]. But in 

that work it was assumed from the outset that / Na represents the probability that the 

J A 

system be in a state characterized by a point of phase space belonging to a domain A in 
a 7n-manifold S'^"' which extends to infinity. Then, invoking that the number of particles 
is conserved, it was concluded that d{Na) is necessarily zero. 

We shall look at these matters in a different way. Let us consider a fixed number n of 
particles, and suppose that their dynamics is given; for instance we write some equations 
of motion. Now we cannot just assume that Na is the above mentioned probability, 
without first proving that its value does not depend on the cutting surface S'^" if the 
same orbits are considered. In other words the definition of "numerical flux of particles" 
should be consistent with the dynamics: its conservation should stem from the equations of 
motion (but we cannot hope to satisfy this requirement without specifying some technical 
information about the nature of these equations). 

The main goal of this article consists in showing that direct interactions, especially when 
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formulated as in predictive mechanics, exactly satisfy this requirement. 

Owing to the advantage of dealing with a finite number of degrees of freedom, a rather 
rigorous setting can be elaborated. We first recollect in Section 2 some useful results of 
differential geometry concerning measure-preserving vector fields which have vanishing Lie 
brackets among themselves (in these matters we try to speak of incompressible vector fields 
only, postponing as far as possible the use of Hamiltonian vector fields). 
Section 3 is devoted to the treatment of direct interactions in the predictive formalism. 
The basic equations of statistical mechanics are introduced following Hakim's early work 
[2]. Then conservation of phase averages easily comes out by straightforward application 
of the geometrical results gathered in Section 2. Not only the fliix of particles but also the 
entropy flux is a closed differential form. 

The Hamiltonian formalism is considered in Section 4; a possible ambiguity related with 
the noncanonical nature of physical positions is discussed and removed with the help of a 
certain idea of the underlying quantum mechanics. Energy, linear and angular momenta 
are introduced therein. 

A possible definition of equilibrium ensembles is tentatively proposed in Section 5, and the 
last Section contains a few concluding remarks. 

2. Geometrical considerations. 

In this section we consider an m-dimensional orientable manifold V"^. Capital indices 
A,B = 1, ....m. 

According to the standard notations of differential geometry, ix denotes the inner product 
(of a differential form) by the vector-field X, and Cx is the Lie derivative. 
In arbitrary coordinates y^, ■■■V^, a volume element is a scalar m-form r/ = 7(y) d"^y, 
where 7 is strictly positive and transforms as a scalar density. If £Ai,...a^ is the Levi- 
Civita object, we define the completely skew-symmetric tensor 

riAi....Am = 7 Am 

Naturally we have set 

d'^y - ^ eA,....Am dy^^ A ....dy^- = dy' A dy^ 

We say that X preserves 77 or equivalently that X is incompressible when CxV = 0. 

For n such that 1 < n < m consider independent vector fields Xi, We associate to 

them the (m — n)-form 

a(X„,,)=^^-^n" d^A,....A^ (2.1) 
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where we define 



dT,Ai....Ar, = 



1 



1lAi,...Ar,Ar,+l...Am 



+1 A ... A dy 



.A' 



(2.2) 



m 



(m — n)\ 



Actually dY^Ax....An is the element of (m — ?i)-surface (imbedded in 1/^"^) induced by r\. 
Indices a,b = 1, ...n. Our notation for a means to recall that a depends not only on the 
vector fields Xa, but also on the volume form. 

Remark: If all the vector fields Xa are tangent to a (m — n)-manifold S'^~'^, then a 
vanishes on this manifold. 

For typographical convenience we define ia = ix^ i = ■ 

In this notation a is obviously proportional to iii2----in rj, which motivates our attention 
to the statements presented below. 

For applications to statistical mechanics it is of interest that a be a closed differential form 
(that is da = 0). Actually this property is ensured by a couple of results presented below. 

Proposition I 

For n < m , consider on the vector Gelds Xi, X2---Xn , preserving the volume element 
rj. If [Xa, Xb] = then for all r < n and for all b we have 



The vanishing of Lie brackets implies [15]that Ca commutes with %. Hence Cahv = 
H^aV = 0. Thus >Cft commutes with all the contractions zi....^- We end up with 
ii ir CbTj, obviously zero. 

We can also assert the following statement 
Proposition II 

Under the same assumptions as above, for all r < n, ii ir r] is a closed form. 

The proof is by induction. diaT] is zero because it coincides with CaV- If the statement is 
true for r — 1, write the identity 



In the second term of the r.h.s. arises the differential of 12... 1^1], which involves r — 1 inner 
products and is a closed form because the statement applies to r — 1. Only the first term 
of the r.h.s. might be nonzero. One is left with 



Ch{ii....ir v) = 



(2.3) 



Proof 



d ii....ir T] = Ci{i2....ir'n) — i\ d{i2....ir v) 



d{ii irv) = C\{i2 irTfi) 
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but this term also is zero according to Proposition I. 

Since a is proportional to ii...i„?7, Proposition II entails that da vanishes as announced. 

Remark The reader who is not familiar with exterior differential calculus can get convinced 
that da vanishes with the help of local coordinates such that Xa = d/dy°' (notice that in 
general 7 7^ 1 in these coordinates). Obviously = 5^, thus = 5^^, etc. Inserting 
these values into a one gets 

« = ^i2....nA„+i....^^ A .... A dy^- (2.4) 

valid only in these particular coordinates. 

With the purpose of applications to averages of various constants of motion, we make this 

statement 

Proposition III 

Under the same assumptions, consider a function f invariant under all the vector fields X, 
that is Caf = Va. We claim that d{fa) = 0. 

In local coordinates where Xa = d/dy"- it is easy to check that df A a — 0. But a global 
statement requires the following calculations. 

Owing to the well-known identity ix{A A B) = ixA A B + {—)'^A A ixB we notice this 
Lemma 

If iaB vanishes for all a — 1, 2....n, then Vr < n, 

{ii....irA) AB = ii....ir{A A B) (2.5) 

easily proved by induction. 

Apply it to ^ = ?7, B = df where Caf = 0. Since degree of 77 is maximal, 77 A = 0, 
which proves that df Aa = 0, hence d{fa) = 0. 

This result will receive several applications in the context of relativistic statistical me- 
chanics: normalization of the distribution function, intrinsic definition of entropy flux, of 
average (linear or angular) momentum. 

2.1. Rescaling 

In connection with the problem of reparametrization of world lines, another property of 

differential geometry is worth noticing. 

Definition. 
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Let Xi, ....Xn have mutually vanishing Lie brackets. An abelian rescaling is the replace- 
ment of these vector helds by 



n 



where (pi, ...tp^ are strictly positive scalars satisfying the condition Ca4>h = for a^h. 
We can obviously check that the new vector helds have mutually vanishing Lie brackets 
among themselves, which justifies our terminology. More general rescalings of the X^s 
would spoil the important Lie bracket condition and will not be considered in this paper. 

Proposition IV 

Let rj be a volume invariant by the vector fields Xi, ...Xn with mutually vanishing Lie 
brackets. Then (0i....0„)~^?7 is a volume invariant by the rescaled vector fields Yi, ...Yn. 
The proof is straightforward: Applied to any form of maximal degree, like the volume 
element, Ca is given by dia. Notice that 



and differentiate. 

Now we immediately get 

Proposition V 

Under the assumptions of Proposition IV, the (m — n)-form a associated with the vector 
fields Xa and the volume element rj is invariant under the replacement X^ (pa^a-, V ~^ 

{(pl....(pn)^^ V- 

3. Statistical mechanics of predictive n-body systems. 

A predictive system of n particles, labelled by subscripts a,b,... running from 1 to n, can 
be defined by the second order differential equations [7] 



This picture exhibits some analogy with the Feynman- Wheeler equations of motion; in 
particular each world-line is parametrized by its own parameter. But in contradistinction, 
the above equations are simply differential equations (they do not involve retarded argu- 
ments). In contrast to Galilean mechanics, the r.h.s. of the equations of motion cannot be 
choosen arbitrarily. This is the price paid for a manifestly covariant formulation. 
The above system is integrable and provides world- lines Xa = Xa{Ta) provided the right- 
hand sides satisfy the predictivity conditions. This scheme amounts to consider the tangent 





(3.1) 
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bundle (T(M'*))"^ as phase space. Setting Va — -j-^ we have in the bundle "natural 

dTa 

coordinates" x,v. The generators of the evolution group are the vector fields 

They can be viewed as linear homogeneous differential operators acting on phase functions; 
in this sense they are Liouville operators. (Accordingly we shall often denote Xf = jCxf 
the Lie derivative of a function / by a vector field X). The lift of world- lines in phase 
space is formed by the integral curves of Xi, X2, ...Xn- The predictivity condition can be 
written in terms of Lie brackets [7] [8] 

[Xa,Xb] = (3.3) 



Notice that these conditions are nonlinear in the "accelerations" ^a- 

The evolution group has n parameters and is abelian; its orbits are the integral curves of 
Xi, ...Xn, and their projection over configuration space M^" gives the world lines. 
A Grst integral is of course a quantity f{x,v) satisfying the relations Xaf = 0,Va. A 
partial integral relative to Xa is defined by X^f = 0,yb^ a. 

When the evolution parameters Ta are affine parameters (the proper-times respectively 
divided by the masses), the r.h.s. of (3.1) satisfy the additionnal condition Va ■ = 0. In 
this case we say that our formulation is autochronous. When the above condition fails to 
be satisfied, we say that the formulation is heterochronous. In this case v'^ are not anymore 
constants of the motion, but they are still required to remain strictly positive. 
Some care is needed however, not only in view of possible ambiguities, discussed later on, 
but also in order to maintain timelikeness of world lines. As soon as condition v ■ ^ is 
relaxed, it becomes necessary to make sure that n first integrals, say Ki, ....K^ can be 
identified as half-squared masses. Assuming that this is actually possible, the physically 
relevant part of phase space is that piece of (T(M^)") which corresponds to positive values 

of Ka. 

Then we can always rescale the Liouville operators to an autochronous formulation, as 
follows. Take 0^ = ^/2Ka/v^. Easy check that >Ca0(, vanishes for a ^ b, which makes 0^ 
admissible for a rescaling according to section 2. If Ya = (j)aXa, and = (f>aVa, it turns 
out that Wa ■ Wa = hence the square of all Wa is constant. In other words, the rescaled 
formulation is autochronous and we end up with n mass-shell constraints = rn^. 
Remark 



8 



Since the Poincare group is implemented in phase space in terms of not only positions 
but also velocities, it is relevant to realize that the definition of this group is invariant 
by abelian rescaling, provided that the scaling factors 0a are themselves scalar Poincare 
invariants. 

Recall this definition : Scalar invariants of the Poincare group [16] are arbitrary combina- 
tions of all the scalar products made of "vector invariants" like: 

Xa - Xb, Vc, {Xa - Xb) AVcAVd 

Why did we bother with consideration of nonaffine parameters? At first sight it seems 
natural to impose affine parameters from the outset, and our first version of predictive 
systems [7]was restricted to this case. But soon it turned out that allowing for arbitrary 
parametrizations has several advantages [17] [18]. Indeed the use of nonaffine parameters 
facilitates the construction of hamiltonian predictive systems in closed form, along the line 
of the a priori hamiltonian approach [18]. This procedure was proved [19]to be equivalent 
with the constraint approach [20]; let us stress that the key for a contact of predictive 
mechanics with this alternative formulation (based upon Dirac's constraints theory) is an 
equal-time condition that is in general not compatible with the affine parametrization of 
world lines. In particular the most tractable toy model of relativistic interaction is a two- 
body harmonic oscillator formulated in terms of parameters which are not affine [8] [18]. 

The structure determined by the vector fields Xi, over phase space allows to formally 

define the evolution group through Uri^____rn = exp (— riXi — TnXn). 

There is a priori no Liouville theorem at our disposal, but we can still look for a conserved 
volume element, that is an everywhere positive Sn-form invariant by action of all the vector 
fields Xa- If such a volume preserved by the motion, say rj, can be found (that is CaTj = 0) 
it is generally not unique. For a given system of vector fields Xa, the invariant volume (if it 
globally exists) can be determined only up to a factor A, which is necessarily a positive first 
integral. Let us write r] = 'j d^'^y, where 7 is some scalar density in arbitrary coordinates. 
Now it is possible to interprete t/ as a unitary operator because we have a scalar product 
for phase space functions, say < f,g >= / f*g rj and this product is invariant under 
evolution. 

It is not yet necessary to assume a Hamiltonian formalism. In fact the contents of Section 
2 imply that several general results can be derived under the simple assumption of a 
conserved volume element; these results seem to open a way to abstract developments in 
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the spirit of mathematical ergodic theory. We shall see later what are the limitations of 
this point of view. 

Asymptotic considerations 

The interest of referring a system of interacting particles to its interaction-free limit has 
been recently emphasized [21]. As temporal conditions have been invoked in order to 
determine this procedure, let us sketch how we can figure the asymptotic time behaviour 
of these particles. It is convenient to remember that a classical scattering theory parallel to 
the quantum one is possible, and has been actually developed in the newtonian context [22] . 
A classical relativistic scattering theory could be considered as well, in close analogy to the 
many-time formulation of quantum relativistic scattering theory available in the litterature 

[23] [24]. Exponentiating the Liouville-operators we get U = exp (— tiXi — t„X„). In 

a similar way the system of n free particles is characterized by the Liouville-operators 

x[^\ Xn^ where = Va ■ da- They correspond to the evolution operator U^^^ = 

exp {—TiX[^^ — TnXji^). Then introduce classical relativistic Moeller operators as the 

limits 

J]± =limt/-^t/^°) 

for Ti, ...Tn alltogether oo. If the limits actually exist in the strong sense, which requires 
that they are independent of the order (as happens in the quantum theory as proved by 
Horwitz and Rohrlich [24]), then it is reasonable to expect the intertwining relations 

un^ = Xa^^ = (3.4) 

This gives the principle of a map between the orbits of the interacting system and the 
solutions of the free-particle motion, invertible map insofar as no bound state is present. 
This seems to suport the claim by Ben Ya'acov that the evolution of the distribution 
function can be fully described in terms of the interaction-free limit of the system [21]. 
However one must be cautioned that no rigorous result about Q,^ is by now available, and 
if anything is going to be proved, it will probably concern short-range interactions rather 
than electromagnetism. 

The Liouville " theorem" is essential * in order to be able of regarding the evolution 
operator U as unitary, in the Hilbert space L^{R^^, rf). In other words one could perhaps 
go without a Hamiltonian structure, but not without a preserved volume element. In 

* Quotation marks refer to the fact that, for the moment, volume conservation is not 
derived from a canonical formalism. 



10 



practice we shall resort to a symplectic form anyway, in order to get rid of the arbitrariness 
of the volume form. 

3.2 Density and distribution function 

Since velocities must point toward the future, the relevant part of phase space is the 
region of T(M^'^) defined hy v ■ v > 0, f'' > 0. The "proper-time-dependent" density 
D{xi, vi, ....Xn, Vn,Ti, ....Tn) is ruled by not one but n Liouville equations 

{^^+Xa)D = (3.5) 

This system of equations was first written by Hakim for a gas of free particles. R.Lapiedra 
and E.Santos pointed out that eq (3.3) just ensures the integrability conditions for n 
Liouville equations in the presence of predictive interactions [11] (unfortunately their for- 
mulation is not thoroughly covariant, which leads them to restrict their statement to the 
first order in the coupling constants). 

As emphasized in ref. [2], only the " proper-time- independent" density 

/+00 
D(xi,Vi, Xn,Vn,Ti,...Tn) dTi....dTn 
-oo 

has a direct physical meaning. ^From (3.5) it follows that AT is a first integral, that is 
XaN = 0. Conversely, eqs (3.5) are formally solved hy D = Ur-^ t„-^5 provided CaN = 0. 

When the masses are specified from the start, the supports of D and are a priori 
restricted to the region 2Ka = m'^; owing to the condition that the velocities point to the 
future, one is left with a one-sheet mass shell. 
Normaliza tion and averages 

The next step will be a manifestly covariant formulation for phase averages and for the 
normalization of the distribution function. We naturally use Hakim's definition [2], that 
is in our notation J Na = 1 where integration is to be performed over a 7n dimensional 
manifold cutting all the orbits of the evolution group. Assuming that the motion generated 
by the Xa 's preserves some volume element rj, we are in a position to apply the results of 
section 2 in the dimension m = 8n. 

Phase space is V^^ — {T(M'^))^ and Xi, ...Xn are the generators of the motion. Now 
m — n = 7n, and rj induces the 7n form 

d^A,...A„ = VA„...A^A^+,...Asr. dy^-+' A ... A 

as element of 7n-dimensional hypersurface imbedded in F^". 
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Since the Liouville operators are incompressible vector fields, we can apply Proposition III 
of Section 2 and conclude that d{fa) = for any first integral /. By the Stokes theorem, 
J fa over any closed 7n-dimensional manifold vanishes. Notice that in particular / can 
be the distribution function. 

We claim that the integral / fa over a 7n-dimensional surface 5^" does not depend on the 
choice of it provided it is specified that S'^'^ cuts once and only once each orbit manifold 
(this precaution discards the possibility of integrating over the mass shell). Since we are 
dealing with a multitime flow, this point is not so obvious as in the one-body case; the 
subsequent details are in order. 

Definition A Row tube is an m-dimensional invariant domain of V^'^ (invariant under 
the n-parameter group with infinitesimal grenerators Xa)- 
Let T be a flow tube. Our precise statement is that 

The integral of fa over the cross-section of T by a In-dimensional surface S'^'^ is inde- 
pendent of 5'^", provided we assume that S'^" cuts (and only once) each integral curve of 

Xi....Xn. 

This can be proved as follows. 

Let T be the orbit of a compact domain A C -S^". The cut of T by another surface 
S"^" will be noted as A'. One must take some care of this complication (absent in the 
one-parameter case): The fronteer dT of the tube has 8n — 1 dimensions, whereas we need 
a 7n-dimensional surface in order to extend the set A U A' as to form a closed manifold *. 
Fortunately it is clear that all the vector fields Xa are tangent to dT. Thus a vanishes 
on dT. Now it is possible to construct a surface B^" imbedded in dT as a submanifold, 
and connecting the domains A and A' in such a way that A U A' U B^"' be a closed 7n- 
surface (see Appendix A) . Then we observe that fa vanishes on B^" for it vanishes on dT. 

Therefore applying the Stokes formula yields / f a = f a. 

J A' J A 

Remark: 

In applications we shall suppose that S"^" = Sf x x 5"^ x (W^)^ where is the space 

of four-vectors and is a spacelike three-surface imbedded in Minkowski space. 
Besides the normalization of A^, another application concerns the total entropy S = 
— J NlogNa. Since the phase function A log A is a constant of the motion, under similar 
assumptions A log Act is a closed form, a result which renders J A log Na independent of 
the integration surface. 

More generally, we have an intrinsic covariant definition for the average of any constant of 
the motion. 

* This complication was pointed out in ref [2] . 
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In general a form like rj is by no means unique. However it may happen, for a given system 
of n interacting particles, that the 8n-form ^d^'^a; A d'^'^v is "accidentally" conserved by 
the motion, providing obviously a preferred choice (see Appendix B for an example at first 
order in the coupling constants). 

Otherwize there remain the possibility of redefining the distribution function (without 
changing parametrization) whenever we change the preserved volume form. So doing we 
keep Na, and therefore the formulas for averages, unchanged. 

Most simply, we shall resort to the Hamiltonian formalism (next section). Although a 
similar ambiguity arises in the Hamiltonian framework, it can be removed by an argument 
which is plausible for all gases made of microscopic "molecules". 

Reduction to 6n dimensions. 

For free particles it was already pointed out by Hakim that normalizing over a 7n-dimensional 
surface actually reduces to a customary normalization over a 6n-dimensional surface, by 
help of the n mass-shell constraints =v'^. It is fortunate that the multitime formalism 
of relativistic dynamics still provides n mass-shell constraints in the presence of interac- 
tions. These constraints amount to fix n constants of the motion Ki, ....Kn identified as 

Wa- 
in most cases of interest, the relevant distribution functions will be concentrated on a 

submanifold 2Ka — m^, with specified values of mi, m^. Therefore we expect that the 

distribution function (as defined in the whole phase space) involves a factor Y[ '^("^a^^-ftTa). 

Since the velocities must point to the future, there is in A?" a factor n^('^a)) where 9 

is the step function. These considerations entail that the effective phase space may be 

Gn-dimensional, as expected for a contact with the newtonian limit. 

Remark: 

Restriction to a sharp mass shell agrees with the realistic picture of a gas of particles with 
given masses; but it is not imposed by our formalism. Moreover, a statistics of particles 
with unspecified masses may be of interest for cosmological applications and for liadronic 
matter. In addition a smearing of the mass shell becomes technically necessary when 
entropy is defined as above, for the logarithm of a Dirac distribution has no mathematical 
meaning. 

In fact the essential mass-shell property is that the quantities — 2Ka are constants of 
the motion. But strictly speaking they are not constraints untill one gives numerical values 
to these constants, as happens when we assume the factors S{rn^ — 2Ka). 

We shall see in Section 4 that the computational complications introduced by the factors 
considered above disappear in the special case of " perfect interactions" . 
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We could hardly go much further without a Hamiltonian formalism; indeed: as soon as 
one whishes to consider either energy or total linear momentum, a symplectic canonical 
formulation is badly needed. 

4. Hamiltonian formalism 

The hamiltonian formulation consists in looking for a symplectic form Q invariant under 
the evolution group, that is Cxa^ = for all a. When an invariant symplectic form Q 
has been determined, the Liouville operators are generated by positive scalars functions 
Ha satisfying = dHa- We refer to these scalars as the Hamiltonians, although they 
are associated to (half) squared masses instead of energy [8]. In any system of canonical 
coordinates QaiPhi we can write Vt = ^dq^ A dpaa- 

The volume form is 77 = d'^'^q A d'^'^p. Notice that r] = const, f]^"^ (exterior power). Since 

= 0, the Liouville theorem follows CaV = 0. 
The scalar density 7 is equal to unity in canonical coordinates; this is not always the case 
in natural coordinates x, v. In arbitrary coordinates, the volume element can be written 

as ri = J d^^v where 7 = ^ ,' , , with an obvious notation for the Jacobian. Drastic 
' ^ ^ ^ D{y) 

simplification arise when the position formulas "accidentally" imply that d^^x d^^p = 

In contrast to Newtonian mechanics, and in view of a famous theorem [25], it is not 
possible to require that the physical positions be canonical variables. In other words we 
cannot have just qa = Xa throughout phase space. This situation implies that, for a given 
system of Liouville operators Xa, the invariant symplectic form Q cannot be unique. This 
complication is a source of ambiguities not only about the Hamiltonian formulation but 
also in the definition of linear momentum (but, insofar as one is concerned only with finding 
first integrals and with solving equations of motion, there is in principle nothing wrong 
with the fact that a given physical system admits infinitely many symplectic formulations) . 
We notice that quantization is naturally affected by this ambiguity, but this is not at all a 
difficulty because there is no reason why a classical system should correspond to a unique 
quantum system; we expect unicity when going from quantum to classical theory and not 
the reverse! 

Remaining in the classical area of physics, the situation becomes more serious as soon as 
we try to construct statistical mechanics with concrete applications in mind. 
If statistical mechanics were limited to formal manipulations involving some preserved 
volume element of phase space, again there would be no problem in having, for a given 
physical system of many particles, infinitely many possible descriptions. But we are inter- 
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ested in the construction of macroscopic quantities like the thermodynamical functions. At 
this stage the Hamiltonian formulation seems to play a crucial role. Several basic concepts 
like the energy or the canonical linear momentum are intimately related with this formu- 
lation through Noether's theorem, and it seems that the definition of most macroscopic 
quantities necessarily involves these concepts; temperature is the first example. Finally the 
very notion of equilibrium rests on the previous determination of a canonical formalism *. 

The lack of unicity in the symplectic formulation can be expressed in saying that the 
dynamics of a rclativistic system cannot be completely determined by the simple knowlege 
of the Hamiltonians. 

However it is true that, to a large extent, spectral and asymptotic properties of relativistic 
dynamical systems are fully determined by the Hamiltonian generators of the motion. For 
example we have been able to speak of the " abstract integration" of a system [8] . Abstract 
integration refers to the solving of the canonical equations of motion, in terms of q,p, and 
provides a foliation of phase space by n-dimcnsional leaves, disregarding the way how each 
leaf is identified as a cartesian products of n lines which are the lifts of world lines. In 
the same spirit, but in the language of constraints theory, it was observed by Todorov 
[20]that, under very general assumptions, the scattering properties of an n-body system 
depend only on the form of the mass-shell constraints, and not on the gauge fixations. 
In other words we can see a relativistic hamiltonian system as an abstract structure (cor- 
responding to the Ha only) completed by additional formulas connecting the physical po- 
sitions with the canonical variables. These position formulas must satisfy the position 
equations [8] 

{x2,Hh}^0 ^a^h (4.1) 
necessary for achieving the determination of world lines. 

To summarize, the ambiguities of the Hamiltonian formalism are of two kinds. 

1) Any change of parametrization implies a redefinition of the "velocities" and a redefinition 
of the set of vector fields Xa- 

2) Even when the parametrization is kept fixed, it is in general impossible to find a unique 
symplectic form invariant by these vector fields. 

* One might speculate that there exists some unifying formula defining temperature, 
entropy and equilibrium, in a way which is invariant under reparametrization and changes 
of symplectic form. But for the moment nobody has any idea of such a formula, if it exists 
at all; therefore one has to cope with the ambiguities. 
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These ambiguities are the price paid for a covariant setting. They cannot arise in Newto- 
nian mechanics for two reasons: the use of an absolute time and the imphcit prescription 
that the physical positions are canonical variables. But in relativistic dynamics, consid- 
ering mutual interactions forbids to require that the physical positions Xa be canonical 
variable. 

At this stage we are led to the following observation. Statistical mechanics is a very general 
theory; its applications may concern a fluid of galaxies as well as an ordinary gas made of 
diatomic molecules. We suggest to distinguish clearly these two situations. For example 
the concept of temperature has certainly a profound physical significance when it concerns 
the air we breath. When applied to a gas of which the "molecules" are galaxies, this 
concept has hardly the same physical meaning, although it is defined through the same 
mathematical structure. Therefore we can provisionnally consider that some ambiguity 
about the statistical mechanics of a very large number of macroscopic bodies is perhaps 
not so scandalous after all. Naturally further work is needed in order to select a convincing 
prescription removing the ambiguity in this case. 

In contradistinction the description of a gas made of microscopic molecules should not 
suffer from the same uncertainty. 

Fortunately, at the scale of an ordinary gas it is reasonable to remember that a classical 
system of n microscopic particles must be considered as the classical limit of a quantum 
system. Our point is that if we carefully formulate the axioms of n-body quantum mechan- 
ics, then taking the classical limit of a system will automatically determine a particular 
Hamiltonian formulation. In order to prove this, let us first sketch the framework for 
quantum mechanics of n interacting particles. The wave equations involve n half-squared- 

mass operators (Hi) op, {Hn)op acting in some Hilbert space and commuting among 

themselves [26] (this point of view stems from quantization of predictive mechanics, but 
it also agrees with constraint theory). However the knowledge of these operators is not 
sufficient for a complete determination of the system. This point is widely overlooked in 
the litterature, because most investigations are concerned with either spectral or scatter- 
ing properties fully contained in the squared-mass operators, all question about position 
measurement being systematically ignored. 

Our proposal was that there exist n coordinate operators with four components, say {x2)op 
satisfying the commutation relations 

[«)op, mop] =0 a^b (4.2) 
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For instance these {x)op^s are relevant objects if one whishes to consider position measure- 
ments *. 

If the system has a classical limit, these relations reduce for ^ — > to the "position 
equations" (4.1) of predictive mechanics [8], and the coordinate operators expected 
to have a unique and well-defined limit x'^{q,p). ^From this limit we can derive the 
generalized Legendre transformation q,p < — > x,v, which also determines the rescaling 
from the parameters ri, ....Tn to the proper times. 

To summarize, the classical limit is expected to determine a predictive dynamical system 
endowed with a preferred invariant symplectic form (corresponding to a preferred Hamil- 
tonian formulation). 

This argument solves an important question of principle. But in practice, the "good" 
symplectic form will remain ignored as long as the details of the underlying quantum 
theory remain ignored. This might be the case of "predictive electrodynamics" [27]untill 
it will be rederived from some n-particle generalization of the Bethe-Salpeter equations. 

4.1. Perfect interactions 

Ideal gases are made of noninteracting particles. The next step consists in considering an 
interaction which can be trivialized by a suitable transformation. With this idea in mind 
let us characterize a perfect interaction by the global existence of particular canonical 
coordinates, say q,p, allowing to write 

Ha = \pI (4.3) 

We refer to q, p as Hamilton- Jacobi (HJ) coordinates and we shall say that we have an 
almost ideal gas. Let us stress the following point: whereas (4.3) could be locally written for 
any Hamiltonian system by solving an n-body relativistic generalization of the Hamilton- 
Jacobi equation, the existence of global coordinates like q,p is a very particular property 
of the gas we consider [28] . 

A system undergoing perfect interactions is spectrally trivial; the Lie algebra of first inte- 
grals has the same structure as in the free system. In principle equations of motion can 
be exactly solved. Nevertheless, word-lines actually deviate from straight lines. It is clear 
that perfect interactions realize the most simple situation beyond the case of noninteracting 
particles. 

* These "coordinate operators cannot, as they stand, be considered as position operators 
in the usual sense. They act on some subset of L^(i?^"', d'^'^x), which is not the case of the 
Newton- Wigner operator. 
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We strongly suspect that considering HJ coordinates amounts to neglect bound states. This 
should not prevent us from using perfect interactions as a simplification when attractive 
forces are present. This approximation may be applied to the case of attractive forces, 
provided that the gas we consider is so dilute that it is reasonable to neglect bound 
states. Naturally one must keep in mind that this representation is valid only insofar as 
the distances between particles are large enough. 

The main interest of H J coordinates lies in the possibility of explicitly writting down covari- 
ant many-body interactions. Indeed, for an almost ideal gas, the most general interaction 
can be explicitly written in terms of the HJ coordinates. 

According to (4.1), the components of Xa must be invariant by all the vector fields 
generated by hamiltonians where b ^ a. We now have to solve {xa-,pl} — 0- 
immediately find 

x: = q:+F^{....p,, q,Apb...) (4.4) 

with all b ^ a m F, and F are arbitrary vector- valued functions except for the following 
restrictions: 

i) They might be required to satisfy reasonable boundary conditions, namely that in 
some sense q,p x,v ai temporal or spatial infinity. 

ii) They must respect Poincare invariance. Scalars invariants and vector invariants of 
the Poincare group can be characterized as well in terms of HJ coordinates. According to 
this remark, in the position formula above must be a vector invariant. 

Of course, the simplest solution x^ — would be trivial; it corresponds to straight world 
lines. 

It might be interesting to discuss under which conditions the Feynman- Wheeler electrody- 
namics can be (approximately) treated as a perfect interaction. 

Using HJ coordinates we can write Xa = p^ d/dq^. For the special choice of a surface S'^^ 
defined by equations of the form 

s{qi) = s{q2) = s{qn) = 

(for instance s{q) = q^) we have this simplification for the 7n-form a of Sections 2,3 

a = const, p^' dcj^, (1) A A da^^ (n) A d^^p + 0(5^^) (4.5) 

where 0{S'^'^) is a term vanishing on 5^" and this notation 

d(J^i{a) = ^£^upa dq!^Adq^A dq^ 
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Then with help of the identity (i(p^) Ad'^p = 2p^ d'^p we factorize d{p^) in d'^^p and finaUy, 
for aU "phase function" $ we can write 

/ = const, i ^ d(pl) A .... Ad(pl) AX (4.6) 

where A(x,ry) is a 6n-form which also depends on the choice made for S"^^. When $ has 

the particular form $ = 9{pfi)....9{p^)5{p^ — m^) 5{p^ — m^)$, then equation (4.6) is 

reduced to 



$ a = const. / $ A (4.7) 

where S^'^ = S'^'^ fl mass shell. This formula makes the contact with more popular 
formulations using a 6n-dimensional phase space; also it is an illustration of a remark 
made in section 3.2. 

It is clear that all first integrals of the interacting system are known in terms of HJ coor- 
dinates, which allows for formally solving the Liouville equations (3.5). But the solutions 
obtained in this way are expressed in terms of canonical coordinates. For practical purpose 
it is usually of interest to have solutions in terms of the natural coordinates x, v, which in 
turn requires knowledge of the position formulas. 

4.2. Average momentum 

In the framework of Hamiltonian formalism, in the absence of external forces the system 
is invariant under translations and rotations, hence the first integrals P = "^Pa, M = 
^ Pa- Application of Proposition III yields an intrinsic definition of the averages 
F , M . hoT mstance 

is conserved under a change of the 7n— dimensional integration manifold. 
5. Canonical Distribution. Equilibrium. 

Provided the system is translation invariant, it is natural to define equilibrium by the 
distribution function 

N = const. 6t 6+ exp{-(3''P^) (5.1) 

for some constant timelike four-vector associated with the inverse temperature. The 
factors 

C = 2^(^a) S{ml - 2Ha) (5.2) 
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take into account the mass-shell constraints and the direction of the arrow of time on world 
lines. 

For an ideal gas, equation (5.1) (which is also an n-body generalization of Jiittner's formula) 
agrees with the definition proposed by Hakim three decades ago *. 

In the presence of interactions, v° is not constant, but 9{v^) is a (discrete) first integral 
although v° is not. Finally AT is a first integral as it should. 

Let us stress that there is in general no indication that the replacement of 0{v^) by 9{p^) be 
legitimate. Strictly speaking, in order to write down the equilibrium distribution function 
in closed form, one might be obliged to solve the position equations! This is no surprise; 
it is reasonable that some information from world lines be necessary. However, interesting 
and drastical simplifications arise in some special will be briefly discussed below. 

Being spatially homogeneous, the equilibrium distribution function is not normalizable, 
which rises the problem of elaborating a covariant procedure for the thermodynamical 
limit. 

In the presence of an external field, space translation invariance is broken, whereas it often 
happens that the energy l^Pfj, remains conserved (the external potentials applied to the 
system must be stationary). In this case {l^Pfi, Ha} = for some constant timelike unit 
vector I defining a "laboratory frame". Then we define equilibrium by 

N = const. 6t....S+ exp(-/3 PP^) (5.3) 

with P a positive scalar. Again, one is led to consider the thermodynamical limit. But the 
analogy with the Newtonian situation is more transparent here. However, existance of the 
limit is an open problem. For instance it was argued that for long-range interactions, it 
may happen that this limit does not make sense [29]). 

Returning to translation invariant systems, it is clear that the natural coordinates x,v 
allow for a simple expression of the factors (at least in the autochronous case); but 7 
may still differ from unity in these coordinates. 

In contrast, any choice of canonical coordonates provides a trivial expression for P, but 
this simplicity is at once destroyed when taking into account the mass-shell constraints. In 
any frame adapted to (3 we get N — const. {Yl ^^) exp(— /3P^), where P° is the sum of the 

* Formula (5.1) seems to differ from eq (6.15) of [2]by the factors 5+. These factors are 
implicit in ref. [2] where the reduction of phase space by mass shell constraints is assumed 
from the outset. 
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p°'s. In ordinary canonical coordinates p° takes on the form p° = ^/m^ + Po ~ 2Va(q',p). 
Also the factors 5+ are in general model dependent, since the canonical expression of Ha 
in terms of q,p can be complicated. 

But in the special case where HJ-coordinates actually exist, remarkable simplifications 
arise. 

First, the mass-shell constraints take on a very simple form (in fact as simple as for free 
particles). In HJ-coordinates = -y/m^ + p2. 

Second, we can legitimate the replacement of by in the step functions. Indeed, 
H = ^jP is a constant of motion, it is clear that on each orbit p remains timelike. Since 
p reduces to v for a vanishing coupling constant, it is clear that p as well as v is oriented 
toward the future. In other words, p and v are simultaneously in the future light cone. 
Third, P = Y,P- 

Finally the equilibrium distribution takes on the free form when expressed in terms of HJ 
coordinates. 

Microcanonical equilibrium 

Let us consider a gas of interacting particles in the absence of external forces. We suggest 
defining the microcanonical ensemble by a straightforward generalization of eq (A4) of ref 
[2], say 

N = const. SiC^" - P^) 5t 5+ (5.4) 

where C is a constant timelike vector and (5+ is given by (5.2). Again, the above distribution 
is a constant of the motion. 

6. Conclusion 

The theory of direct interactions provides a solid ground for a covariant formulation of 
statistical mechanics. Let us be more specific: predictive relativistic mechanics is the 
most natural way for taking seriously the idea (already present in Hakim's early work) 
that one is dealing with an n-parameter evolution group, associated with a multitime 
parametrization. As we know for a long time, this philosophy can be expressed in terms of 
modern differential geometry, which leads to vector fields in 8n-dimensional phase space; 
these "Liouville operators" have mutually vanishing Lie brackets, ensuring that initial 
velocities and positions actually determine world lines. 

If we now consider a statistical mechanics of relativistic particles, the predictivity condi- 
tions are technically essential; on the one hand they entail integrability of the n Liouville 
equations (it is an exact result, by no means limited to first order in the coupling); on the 
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other hand they ensure, as we have shown, the intrinsic nature of important quantities 
defined by a flux integral. For instance they imply the conservation of phase average val- 
ues, especially the conservation of norm (numerical flux of particles) and entropy. A few 
manipulations of exterior differential algebra were sufficient to demonstrate that predictive 
mechanics exactly fulfills the demands of relativistic statistical mechanics concerning the 
Lorentz invariance of the phase averages. 

In our exposition we have as much as possible separated the specific features based upon 
the Hamiltonian formalism from more general topics requiring simply a preserved volume 
element. 

Insofar as one aims at a general theory devoted to a description of systems out of equi- 
librium, the basic notion is that of a preserved volume element, rather than a symplectic 
structure: in other words the Liouville theorem without the Hamiltonian formalism! This 
viewpoint is close to the spirit of mathematical ergodic theory. However we notice that, for 
a system defined by its equations of motion, the preserved volume element is not univocally 
determined. 

In this work we have focused on the n-particle distribution function. But we expect that 
also reduced distribution functions can be proved to satisfy (in the presence of mutual inter- 
actions) covariant evolution formulas and will enjoy conservation properties with respect 
to normalization. This point is obviously related with the important issue of a tractable 
and relativistic BBGKY hierarchy. 

A Hamiltonian formulation is highly desirable anyway, as soon as one is interested in a truly 
physical description. In the presence of mutual interactions, this formalism is essential in 
order to define linear momentum and energy. Thus the notion of equilibrium and probably 
also the construction of thermodynamical functions require a Hamiltonian formulation. 
At least for gases made of microscopic constituents, a certain idea of the underlying quan- 
tum mechanics solves (in principle) the ambiguity related with the impossibility of having 
canonical physical positions. 

The usual expression for the equilibrium distribution function has been generalized in a 
covariant formula taking mutual interactions into account. Further work is needed in order 
to check in which sense and in the context of which side conditions this definition actually 
realizes the extremum of the entropy integral. 

There is a good hope that several methods and formulas concerning ideal gases remain 
valid when the mutual interactions are modeled as "perfect interactions". This simplified 
picture provides a toy-model of interacting particles for which rigorous results can be 
proved. Its relevance for a real gas is probably limited to the dilute-gas approximation. 
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but many standard calculations of the theory of ideal gases can be extended to almost 
ideal gases in a straightforward fashion. 

Rather than new practical applications, the scheme proposed here opens theoretical pos- 
sibilities. For instance it sheds light on the transformation properties of various results 
available untill now in noncovariant form. 

The present formulation allows for considering particles with unspecified masses, which 
is interesting for applications to cosmology or to hadronic matter (the latter suggests an 
extension toward quantum relativistic statistical mechanics). 

Several generalizations are possible; for example one could enlarge phase space as to include 
radiation. 

As it stands our exposition may seem to be exceedingly geometrical, but this feature is an 
advantage in view of possible generalization to curved spacetime. 

A comparison with the "statistics of events" developed by Horwitz et al. [13]will be of 
interest. 

Appendix A 

Through any y G A passes a unique orbit, say of the evolution group. According to the 
assumptions made, this orbit is cut in a unique point y' by the surface 5"^". We can write 

y' = Ur'^r'2....T'r. y 

for a sequence of parameters r'l, r'2, t' n which depends only on y and on the cutting 

surface. Among all possible curves lying in F^, a very simple way for joining y to y' is 
characterized for instance by the conditions 

Tg ^ r'a 
Tl r'l 

Indeed F is endowed with a network made of the integral curves of Xi, This proce- 

dure defines a distinguished curve Cy which lies in Fj, and joins y with y' . 
Now when y runs in the {In — 1) -dimensional fronteer of the domain A, the curves Cy 
generate a 7n-dimensional surface connecting 5 A with 5A'. Notice that two different 
F cannot intersect unless they are identical. When y runs in 5A the n-dimensional orbits 
F generate dT. Since each Cy lies in Fj^ it is clear that, in the generic case, B^"' C dT as 
a submanifold. 

In this construction, B^"' can be considered as the fronteer of a "sub-tube" generated by the 

integral curves of a vector field S = Xi + V2X2 + + J^n^n where the n — 1 independent 

quantities 1/2, ■■■■J^n are first integral of the system Xi, 
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Appendix B 

At first order in the coupling constant, the 8n-form 77 = d'^'^x A d'^'^v is preserved by- 
predictive electromagnetism and by predictive vector interaction. 
Proof 

At this order of approximation, the acceleration on particle a is given by [27] 

where the skew-symmetric tensor Fa^b-, is (up to a constant factor) the Lienart-Wiechert 
field created at point Xa by the (vectorial or) electric charge with label b supposed to 
move on a straight line in spacetime. Therefore, Fa^b depends on the VbS for b ^ a, but 
certainly not on Va- Hence d^"/dv'^ = F^. By contraction of indices d^^/dv^ = 0. Since 
we are using the natural coordonates x, v, this equation amount to write OaX"^ = 0, which 
proves our statement. 

The author is grateful to Prof.Remi Hakim for a stimulating discussion. 
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